Thermal rectification and negative differential thermal resistance in lattices with mass gradient 
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We study thermal properties of one dimensional(lD) harmonic and anharmonic lattices with mass gradient. 
It is found that the temperature gradient can be built up in the ID harmonic lattice with mass gradient due to 
the existence of gradons. The heat flow is asymmetric in the anharmonic lattices with mass gradient. Moreover, 
in a certain temperature region the negative differential thermal resistance is observed. Possible applications in 
constructing thermal rectifier and thermal transistor by using the graded material are discussed. 
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Theoretical studies of heat conduction in low dimensional 
nonlinear lattices in past years have not only enriched our un- 
derstanding of the microscopic physical mechanism of heat 
conduction, but also suggested some useful thermal device 
models such as rectifier/diode |1] and thermal transistorlQl 
for controlling heat flow. More importantly, a two segment 
thermal rectifier has been realized experimentally by using 
nanotubeslSJ, which implies that we can control and manipu- 
late phonons like we do for electrons. 

On the other hand, the functional graded materials (FGM) 
have attracted increasing attention in many fields ranging from 
aerospace, electronics, optics and bio-materials etc. due to the 
intriguing properties |4, 5]. The FGM are a kind of inhomoge- 
neous materials whose compositions and/or structures change 
gradually in space which results in corresponding changes in 
physical properties such as electric, mechanical, thermal, and 
optical properties. The FGM are abundant in nature, and can 
be also purposely manufactured fll^. 

However, compared with optical, mechanical and many 
other properties, the thermal properties of the graded mate- 
rials have not yet been fully studied (see recent review article 
Inland the references therein). 

In this paper, we study thermal properties of the FGM rep- 
resented by ID harmonic and anharmonic lattices with mass 
gradient. For convenience, we call the lattices mass graded 
harmonic/anharmonic lattices. They can be used to model su- 
perlattice or layered structures. As it will be seen that the Id 
graded lattices exhibit new physics such as asymmetric heat 
flow and negative differential thermal resistance, two essen- 
tial properties used to construct thermal rectifier and thermal 
transistor. Therefore, the graded material might find new ap- 
plications in heat controlling and management. 

We consider a ID mass graded chain, which is equivalent 
to a chain with graded coupling constants|7]. FiglTJa) shows 
the configuration. The mass of the particle is M, = Mmax - 
(i - 1){Mmax - Mmin)/{N - 1), where Mmax is the mass of 
the particle at the left end and Mmin is that of the particle 
at the right end. is the total number of the particles. The 
Hamiltonian of this model is: 



here jc, is the position of the particle. Without loss of gener- 
ality, V takes an anharmonic form, namely, (jc, - - fl)^/2 + 
/3(xi - Xi-i - fl)'*/4, which is a Fermi-Pasta-Ulam (FPU) /3 po- 
tential. So the system is called a graded FPU lattice. The 
lattice constant a = 1 .0 and the coupling constant f3 - 1 .0. In 
the case of y6 = 0, this lattice reduces to a graded harmonic 
chain. 

In our simulations we use both fixed and free boundary 
conditions. The stochastic heat baths are put on the 1", 2"'', 
{N-l )''' and A^''' particles with temperature Tl and Tr respec- 
tively. The equations of motion (EOM) of these four particles 
are described by the Langevin equations: 



MiXi = F(xi-Xi^i)-F(xi+i-Xi)- 



- '^L^i, i — 1,2 

4R-ARXi, i^N-l,N, 
(2) 

where ^l/r are independent Wiener processes with zero mean, 
waiiancelAi/RkBTi/R and force F = -dV/dx. To minimize the 
Kapitza resistance between the bath particle and its neighbor, 
we set Ai/R/Mi = 0.1. Both the fourth-order Runge-Kutta and 
the velocity Verlet algorithm are used to integrate the EOM. 
Differences between results of these two integration meth- 
ods are negligible. Simulations are performed long enough 
(> 10^ time units) such that the system reaches a stationary 
state where the local heat flux is constant along the chain. 

We start with the mass graded harmonic lattice, i.e. y6 = in 
the potential V(x). FiglTfb) shows the temperature profile with 
different lattice length N and different boundary conditions. 
Mmax and Mmin are fixed at 10 and 1, respectively. It is well 
known that no temperature gradient can be built up along a ho- 
mogenous harmonic lattice|8]. However, in the mass graded 
harmonic lattice, temperature gradient is clearly seen. 

The local flux at site / is defined as 



Ji = ^fl<(i,+i + Xi)F(xi+i - Xi)). 



(3) 



H 



(1) 



After the system reaches a stationary state, 7, is independent 
of site position /, so the flux can be denoted as J. The heat 
conductivity is calculated by /c = -J/(dT/dx). We should 
stress that in calculating the temperature gradient, we get rid 
of temperature jumps at the two boundaries. In FiglTJc), the 
heat conductivity versus N for different boundary conditions 
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FIG. 1: (Color on-line) (a) Schematic picture of the mass graded 
harmonic lattice, (b) Temperature profile for the lattice with fixed and 
free boundary conditions, (c) The heat conductivity k versus system 
size N for the lattice with Mmax = 10 and M^in = 1. k = AN, 
where A = 0.048 and 0.027 for free and fixed boundary conditions, 
respectively. The temperatures of two baths are = 1.1 and Tn = 
0.9 for left and right end, respectively. 



are plotted. For the same lattice length A^, heat conductivities 
with fixed boundaries are lower than those with free bound- 
aries, because there is a limitation that all vibrational eigen- 
modes must vanish at boundaries in fixed boundary cases. It 
is clearly seen that the heat conductivity diverges linearly with 
length A^. This linear property is independent of boundary 
condition. It is different from that of disordered harmonic 
lattice where the thermal conductivity diverges with system 
size as and 1 / for free and fixed boundary condition, 
respectively jH. 

The linear divergence of thermal conductivity on system 
size is a new phenomenon. As we know that in ID homoge- 
nous harmonic lattices, no temperature gradient can be built 
up and thus thermal conductivity cannot be defined ! lOll . In 
order to understand the underlying mechanism, we have to in- 
voke the vibrational eigenmodes of the graded system. It is 
found in Ref.|7] that there is a critical frequency coc, which is 
the maximum eigenfrequency of the coiTesponding homoge- 
neous harmonic lattice with M,- - Mmax, where Mmax is the 
maximal mass in the graded harmonic lattice. The modes with 
a> < a>c can be well extended over the whole chain, whereas 
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FIG. 2: (Color on-line) The heat conductivity k versus system size A' 
for graded FPU lattice under free boundary conditions (a) and fixed 
boundary conditions (b) with Mm^x = lOandM^/A, = 1. In all cases, 
K oc N". The values of a are given in the figure. 



those with oj > ojc ( called gradons) are localized at the side 
with lighter masses. Therefore, the formation of temperature 
gradient in the graded harmonic lattice is attributed to the lo- 
calization of the gradons. 

In the following, we focus on the mass graded anharmonic 
lattice. We first examine the size effect of heat conductiv- 
ity and its dependence on temperature. It is found that the 
thermal conductivity diverges with system size as, k ~ N". 
Fig|2la) shows the temperature effect of the divergent expo- 
nent under free boundary condition. The value of a does not 
change very much, it almost keeps the same value of 0.36 
(0.35) when the average temperature To is increased from 1.0 
to infinity. The value of the divergence exponent is very close 
to the results from the renormalization group for Id hydrody- 
namic systems OI III , and that one for the hard core model lll2ll . 
The calculation for infinite To is actually realized by discard- 
ing the quadratic term (of potential) in Hamiltonian since this 
term is negligible compared with the quartic term in infinite 
temperature limit. 

Fig|2|b) shows the temperature effect of the divergent expo- 
nent under fixed boundary conditions. In the case of To - I, 
the best fitting with all available date up to - 512 gives 
rise to an a - 0.51. However, the best fitting for values of 
the largest five A^ (512 < A^ < 8192) gives rise to a = 0.35 
which is very close to the value at the infinite temperature. 
The a - 0.51 might be a finite size effect. 

The value of a for free boundary condition is quite close to 
the value for fixed boundary condition. Therefore, this result 
seems to be very similar to the homogenous FPU-/? lattice, 
namely, the divergent exponent a seems to be independent on 
the boundary conditions. 

As for the disordered FPU-yS lattice, it is observed that the 
value of jS depends on the boundary condition lfTil . How- 
ever, this result needs to be further checked as the authors in 
Ref.|fl3l| used very short lattice, A^ < 20. 

We should point out that the above results may not be very 
conclusive as we have an intermediate value of system size A^. 
In order to get a more convincing conclusion, one needs to go 
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FIG. 3: (Color on-line) (a) Normalized Heat flux Jn versus A for 
three different mass gradients. A' = 200, To = 0.1, M^^™ = 1 and 
Mmax = 20, 10 and 5. Jn = JIC, where C = 0.0045, 0.0080, and 
0.0145 for Mmax = 20, 10 and 5 respectively, (b) Normalized Heat 
flux Jn versus A for To = 1.0,0.1, and 0.01. A' = 200, M„™ = 
^,Mmax = 10. Jn = JIC, where C = 0.072,0.008 and 0.0005 for 
To = 1.0,0.1, and 0.01 respectively. 
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FIG. 4: (Color on-line) (a) Heat conductivity k versus temperature 
To. A = 0.1. (b) Heat flux J versus T,. T, is fixed at 0.4. A' = 
200, Mmin = 1 and Mmax = 20 and 10. The negative differential 
thermal resistances (R < 0) appears when Ti < 0.14 and 0.11 for 
Mmax = 20 and 10 respectively. 



to large N, say > 10^ which is very difficult for current 
computer facilities. However, this is not the main purpose of 
the paper. 

In the following, we study the thermal rectification in the 
mass graded anharmonic lattice. In order to avoid the ef- 
fect of the Kapitza resistance at the boundary, we record the 
temperatures of the closest particles to bath particles at two 
ends, namely, the and the 198"' for A^ = 200 as T,,,. 
For convenience, we set the temperature of heavy mass end 
T, = Tai l + A) and that of light mass end T,- = ro(l - A), 
where Tq is the average temperature of the system. 

For comparison, we show the normalized heat flux Jn{= 
J I C) versus A in Fig|3ja) for three different mass gradients, 
where C is a constant. When A > for different mass gra- 
dients, the heat flux increases steeply with the increase of A. 
However, in the case of A < 0, the heat flux changes a little 
when A changes. This asymmetry of heat flux with respect 
to A is called thermal rectification. We study the dependence 
of the rectification on the mass gradient (fix Mmin — 1 and 
change Mmax - 20, 10 and 5). It is shown in Fig|3la) that the 
larger the mass gradient, the more obvious the rectification. 

To find the temperature dependence of rectification, we cal- 
culate normalized heat flux Jni versus A for different temper- 
ature, Tq = 0.01,0.1 and 1.0, (Figl^b)). From this figure 
we can see that the thermal rectification vanishes in both high 
and low temperature limits. In the low temperature limit, the 
graded anharmonic lattice reduces to the graded harmonic lat- 



tice in which no rectification exists. In the high temperature 
limit, the low frequency vibration modes, which dominate the 
heat conduction, in both ends can be excited, therefore no rec- 
tification is found either In the case of To = 0.1, the low 
temperature end is harmonic, while the high temperature end 
is strongly anharmonic, thus the vibrational spectra strongly 
mismatch with each other, which leads to the thermal rectifi- 
cation. 

Finally, we should stress that the rectification is sensitive 
to boundary conditions. The effect can be observed only 
for fixed boundary conditions, because of the low frequency 
mode cannot be restrained for free boundary conditions. 

We show temperature dependence of heat conductivity in 
FiglUa) for the graded anharmonic lattice. Another important 
feature found in this ID graded chain is the negative differ- 
ential thermal resistance (NDTR)^?], namely, the larger the 
temperature difference the less the heat flux through the sys- 
tem. In order to illustrate the NDTR, we fix the temperature 
of the right (light mass) end at Tr - 0.4 and plot the heat flux 
J versus temperature Ti in FigHJb). The differential thermal 
resistance is defined as /? = ^^^J I ifj^^consi- NDTR is 
seen on the left part of the vertical line, where J increases as 
Ti is increased. 

To understand the mechanism of rectification and NDTR, 
we calculate the power spectrum of particles close to the two 
ends, and then compare their spectrum with each other(Fig|5]l. 
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FIG. 5: (Color on-line) Spectra of the two end particles for system 
A' = 200, To = 0.1, Mmw = 1 and Mmax = 10. (a) The heavy mass 
end (left) contact with high temperature bath (A = 0.6) (b) The light 
mass end (right) contact with high temperature bath (A = -0.7). 



found that the temperature gradient can be built up in the ID 
graded harmonic lattice chain due to the localization of high 
frequency gradon modes. The heat conductivity diverges with 
the system size linearly, k cc N. In the graded anharmonic 
lattices, the thermal conductivity diverges with system size as, 
K ~ N", the value of a seems to be independent of temperature 
and boundary condition. 

The thermal rectification and the negative differential ther- 
mal resistance have been observed in the graded anhar- 
monic lattice. This is quite similar to the recent nanotube 
experimentfs'l, in which half of the tube has been gradually 
mass-loaded on the surface with heavy molecules. Our results 
suggest that the graded materials might be used as thermal 
rectifier and thermal transistor 
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When temperature is high enough, quartic term in FPU poten- 
tial is the dominant term in the whole chain. So the coupling 
among modes is strong and low frequency modes, which con- 
tribute to the heat conduction most, are abundant in the spec- 
trum. Then the flux depends mainly on temperature differ- 
ence. So the quartic term is only dominant at the end with 
high temperature and the quadratic term plays the main role at 
the low temperature end, therefore, the low frequency modes 
cannot go through the system. The heat current is controlled 
by two competing effects: temperature gradients and overlaps 
of vibrational spectra. In the presence of mass gradient, there 
is a big difference between A > and A < in the spectra. As 
shown in FiglSja), where the left end (heavy mass) contacts 
with high temperature bath (A > 0), the spectra of the two 
particles overlap in a wide range of frequencies, thus the heat 
can easily goes through the lattice along with the direction 
of temperature gradient. Whereas when the right end (light 
mass) contacts with high temperature bath (A < 0,Fig|3b)), 
the spectra separate from with other. It can be seen that the 
light mass particle with high temperature oscillates mainly 
in high frequency, however the heavy mass particle with low 
temperature oscillates in low frequency. As a result, heat is 
difficult to go through the system although there is a temper- 
ature gradient. When Tq goes to low temperature limit, the 
graded chain can be regard as a graded harmonic lattice whose 
thermal properties are like those of the homogenous harmonic 
lattice, then the flux depends only on the temperature differ- 
ence. 

In summary, we have studied thermal properties of ID mass 
graded harmonic and anharmonic lattices in this paper. It is 
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